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Abstract 

In this paper we show that a process modeled by a strongly con- 
tinuous real-valued semigroup (that has a space convolution operator 
as infinitesimal generator) cannot satisfy causality. We present and 
analyze a causal model of diffusion that satisfies the semigroup prop- 
erty at a discrete set of time points M := {r m | m € No} and that is in 
contrast to the classical diffusion model not smooth. More precisely, 
if v denotes the concentration of a substance diffusing with constant 
speed, then v is continuous but its time derivative is discontinuous 
at the discrete set M of time points. It is this property of diffusion 
that forbids the classical limit procedure that leads to the noncausal 
diffusion model in Stochastics. Furthermore, we show that diffusion 
with constant speed satisfies an inhomogeneous wave equation with a 
time dependent coefficient. 



1 Introduction 

The standard model of diffusion and its variants have many applications, for 
example, in solid state physics ([II]), environmental modeling (IS]), image 
processing ( jBJ EH E21 E3] ) and inverse problems ( [3] m El 031 EI] ) • For these 
applications the standard diffusion model ([6]) 

dv 

(1) — - V • (D Vv) = t > with v\ t=0 = u 
at 

is usually modified by replacing the diffusion constant D by a tensor that 
takes inhomogenities or anisotropy of the medium into account. Recently sub- 
diffusion (very slow diffusion) and super- diffusion (very fast diffusion) have 
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been modeled by a completely different approach where the time derivative 
in (TjQ) is replaced by a fractional time derivative (see e.g. [12] )• Although, 
many of these models (or their discretizations) are successfully used in ap- 
plications, they violate the principle of causality. By causality we mean, for 
example, that a characteristic feature of a process like an interface or a front 
must propagate with a finite speed. It is peculiar that causality is a funda- 
mental concept in physics of wave propagation and the respective mathemat- 
ical theory of hyperbolic equations, but it is ignored in physics of diffusion 
and the respective mathematical theory of parabolic equations and stochas- 
tic processes. Is it so difficult - or even impossible - to find a causal model 
of diffusion? It seems that processes like diffusion are more "complicated" 
than processes related to waves, but as a matter of fact many of the stan- 
dard models of wave propagation in dissipative media (cf. [ISlEIUdTlITS]) ^ ac ^ 
causality, too. Causality demands a coupling of space and time and therefore 
forbids model equations that are to "simple". Hence, apart from its physical 
significance, it is interesting to study causality from the mathematical point 
of view. Even if causality is not relevant for many applications, which is not 
really clear yet, there is no mathematical reason to ignore this issue. 

It is the concern of this paper to model and analyze a causal model of 
diffusion. The paper consists of essentially two parts. First we show that 
causal diffusion cannot satisfy the strongly continuous semigroup property 
and its respective evolution equation. More precisely, we show that a process 
decribed by the scalar-valued evolution equation 

Ov 

(2) Av = t>0 with v\ t=0 = u, 

at 

where the space convolution operator A is the infinitesimal generator of a 
strongly continuous semigroup (0 [5]), cannot satisfy causality. That is to 
say, there exists a fundamental solution G : ~R N x K — > R of (J2J) that does 
not satisfy 

(3) supp(G) C {(x,t) eR N x [0,oo) | |x| < c t} 

for a positive constant Cq. This causality condition can be interpreted as 
follows: Consider the diffusion of an ink droplet of radius Rq dropped in 
water at time t = 0. Then causality demands that the ink water interface 
propagates with a finite speed bounded by Cq. Even if there is no interface 
(visible), the concentration v of ink must vanish outside a closed ball of radius 
R(t) <R + c 1. 

Although diffusion equation ([1]) does not satisfy ([3]), many scientists con- 
sider this equation as weak causal (or shortly as causal), since there exists a 
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fundamental solution G that satisfies 



(4) G{-,t) = Q for t<0. 

However, this property has nothing to do with causality. As a consequence, 
the problem of causality has beed ignored. 

The second and larger part of this paper presents and analyzes a causal 
model of diffusion. In case of constant speed Co of diffusion, the model reads 
as follows 

f a m b fx', n(t) t) dcr(x') 

\p) V(x.,t) = j— With v\ t=0 = U, 

\ b R(t){U)\ 

where r > denotes a small time period, n(t) G No is such that t G 
(nT,(n + 1) r] and R(t) := Co (t — n(t)r). Moreover, dcr(x') denotes the 
Lebesgue surface measure on WL N and 15^(0)1 denotes the surface area of the 
sphere Sr(0). This diffusion modell satisfies causality, the semigroup prop- 
erty on the discrete set of time points M = {m r | m G No} and i 4 |(x, t) 
is discontinuous at t G M for each x G supp (v(-,t)). Thus model (jHJ) consti- 
tutes a compromise between causality and the strongly continuous semigroup 
property. Physically, this model describes a transport process of mass during 
which the masses split up at time points t G M and then spread out in each 
direction without interfering!]] As a consequence, not each space-time point 
has a unique velocity vector, but the general continuity equation 

dv 

(6) — (x,t) + V-j(x,t) =0 for t^mr, m G N 
at 

holds with flux density 

• < \ f v(x + R(t)y,n(t)r) , , N 

(7) !{*.*) = -C J |s l( 0)i ydg(y) - 

Si(0) 

For those coordinates (x, t) for which 

u(x + y, r) w u(x, r) + i2(t) V«(x, n(t) r) • y 
holds, it follows that 

c Qj R(t) 



(8) j(x,t) « -D(t) Vv(x,n(t)r) with !>(*):= 



lr The superposition of diffusion processes over a continuous range (0, c*] of speeds does 
not satisfy the strongly continuous semigroup property, but the splitting up of masses take 
place at each time point t > r* for some r* > 0. 
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since 




ykViMy) = 8k,i \Si(0)\/N 



k, I e {1,2 



...,N}. 



Si(0) 



But this approximation of the flux density resembles Fick's lavJ^ 



D Vv(x,t) 



with 



D := D(t/2) . 



Note if r — > say under the side condition D (r/2) = const ., then Co — > oo, 
i.e. causality does not hold for this limit. 

Finally, we show that the diffusion process §J5§ satisfies an inhomogeneous 
wave equation with time dependent coefficient such that mass is conserved. 

This paper is organized as follows: In Section [2] we prove under fairly 
reasonable assumptions that causal diffusion cannot satisfy the property of 
a strongly continuous semigroup. For the convenience of the reader we put 
some of the technical parts of the proof of the main theorem in the Appendix. 
Subsequently we present and analyze a causal diffusion model in Section [3J 
The paper concludes with a short section of conclusion (Section HJ). 

2 Real-valued semigroups and causality 

It is common in science to assume that the operator A in (j2J) has nice prop- 
erties like linearity, translation invariance and some kind of smoothness such 
that it can infered that A is a space convolution operator. Theorem 4.2.1 
in [TO] provides such a list of necessary properties. For example, translation 
and rotation invariance of A means that the medium in which the process 
takes place is homogeneous and isotropic. In this section we assume that A 
is a space convolution operator and an (infinitesimal) generator of a strongly 
continuous semigroup in L 1 (IR 7V ) and prove that the solution of evolution 
equation (J2J) cannot satisfy causality condition (j3J) • This result suggests that 
even if some of the nice properties of A are dropped, causality cannot be 
restored. And therefore we present a new model for causal diffusion that do 
not satisfy the semigroup property for all time points in the following section. 

For the convenience of the reader we put the technical parts of the proof 
of the main theorem in the Appendix. 



2 Cf. Remark Hat the end of Section [3 
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2.1 Description of the evolution problem 

The notion of the Fourier transform and the convolution used in this pa- 
per are specified at the beginning of the Appendix. Moreover, we use the 
following notion 

Definition 1. Let N G N and a : R — > R be a positive and rotational 
symmetric function such that 

exp(-a(k)) G S(R N ) and (2n)- N ' 2 J | J- 1 {exp(-a(k))}(x)| dx = 1 . 

Moreover, let a := (2 ti)- n ' 2 ^{a}. The operator A : S'{R N ) -»■ S'(R N ) is 
defined as the convolution operator 

(Au)(x) := -(a* x u)(x) for u G S'(M. N ) . 

Here * x denotes the convolution with respect to the space variable x G M. N . 

This definition is reasonable, since the convolution / * x g of distributions 
/ G S(R N ) and p G «S / (R JV ') is well-defined and lies in S'(R N ) (cf. e.g. [Zl DZ3 ) - 

Let a and .4 be as in Definition [TJ Then 

(10) G(x,i) := (2 7 r)- iV/2 ^- 1 {exp(-at)}(a;) (x G R, t G R) 
is well-defined and defines by 

(11) (St«)(x) = *x«)(ac) foruG5'( 



a semigroup {5 t : 5 / (R JV ) — >■ ^'(R^) 1 i > 0} with generator A (cf. e.g. 
Chapter XVII in pj or Chapter 7.4 in [5]). Moreover, G satisfies the evolution 
equation 

(12) ^_^ G = on fx(0,oo) 
with initial condition 

(13) G(x,0) = 5(x). 

Here <5(x) denotes the dirac delta distribution on R^. Conversely, if {S t \t > 
0} is a semigroup with a generator A that is a space convolution operator 
with kernel a, then G defined as in (flUD satisfies (ITT1) and ffT2l) with ffT3l). 
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Remark 1. The assumption that a is rotational symmetric is motivated by 
the following fact. If a process takes place in a homogeneous and isotropic 
medium, then - by symmetry - the fundamental solution Gq of the respectively 
evolution equation must be rotational symmetric in x. And consequently, a 
is rotational symmetric, too. Since this paper is a first approach on causality 
of semigroups and its evolution equations, we focus on these special types of 
processes. 

Example 1. A prominent semigroup is defined by the Gauss function 

(14) G(x, t) = (4tt D ty n ' 2 exp {~^\ (x G M. N , t G R), 

where D is a constant. The respective generator A is defined by 

o(k) := D |k| 2 (kGM w ). 

The process described by this semigroup is not "causal", since the Gauss func- 
tion G is everywhere positive. Moreover, G satisfies the standard diffusion 
equation (T7]j with w(x) = S(x). □ 



2.2 The structure of a 

We first investigate the properties of a, the Fourier transform of the kernel a 
of the operator A. For this purpose we recall some facts about holomorphic 
roots. 

If / and g are entire functions such that / = g 11 , then we call g an n—th 
holomorphic root of /. If the entire function / has an n— th holomorphic 
root and 

(15) /W(c) = 0, / (1) (c) = 0, f^\c)=0, / (m) (c)^0. 

for some c G C, then n divides m. Here f^ := / and f^ 1 ' denotes the n— th 
derivative of f(z) with respect to z. For a proof of this statement see the 
"Wurzelkriterium" in Chapter 3 Section 1 in [20J. 

Theorem 1. Let f be an entire function. If f has no zeros, then f has 
an n—th holomorphic root for each n G N. Conversely, if f has an n—th 
holomorphic root for each nGN, then f has no zeros. 

Proof. The first statement is well-known (cf. for example the "Wurzelsatz" 
in Chapter 3 Section 2 in [19J). For the second statement. Assume that / 
has a holomorphic root for each n G N. Then for each c G C there exists a 
finite m G N U {0} such that (IT5|) holds. According to the previous remark, 
each n G N divides m and consequently m must be zero. But this means 
that / has no zeros which concludes the proof. □ 
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Now, with the help of Theorem [T] and the Paley- Wiener- Schwartz Theo- 
rem, we can prove that z G t— > a(z) is entire if (J3]) holds. 

Lemma l.IfG defined by ( f77J|) satisfies condition then a can be extended 
to an entire function denoted by z G C i— >■ a(z). 

Proof. Let 

(16) / t (k) := T{G(-,t)}(k) for k G R N , t > . 

If G satisfies ([3]), then according to the Paley- Wiener-Schwartz Theorem 
(cf. Theorem [5] in the Appendix), each /t(k) can be extended to an entire 
function, denoted by F t (z), such that 

(17) F t (k) = f t {k) = exp{-a(k) t) for keR N . 

Moreover, it follows that Fi/ n (z) is an n— th holomorphic root of Fi(z), i.e. 
(Fx/ n (z)) n = Fi(z) and Fy n (z) is entire. This together with Theorem [p, 
implies that F\ has no zeros. Hence there exists an entire function a (z) 
such that 

F t (z) = F t (0) exp{-a (z) t} with a (z) = - J d^ , 

where 7 2 is any path connecting and z. Since (fTTj) holds, it follows that 
a (k) = a(k) for k e R and -F<(0) = 1. This proves the lemma. □ 

2.3 Proof of noncausality 

For the convenience of the reader we recall some notion about polynomials 
in several variables. Let z = (zi, . . . , Zfj) G and n := (nx, . . . , n^) G N^. 
Then 

N 

(18) z n := z^ . . . z n N N and \ z \ 2 :=^z 2 n . 

n=l 

In the following chsupp(g) denotes the closed convex hull of the support of 
9- 

Theorem 2. Let a and A be defined as in Definition^ and G be defined as 
in / TJOj) . Moreover, let f\ := G(-, 1). If a is entire and 

(19) chsupp(fi) = B CQ (0) with Co = 1 , 
then supp(Afi) C B CQ (0) cannot hold. 
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Proof. We perform a proof by contradiction and assume supp(^4/i) C £>i(0). 
There are two cases, either J 7 ' 1 (a) has or has not compact support. 

i) Assume that a has compact support. Then according to the support 
theorem for distributions with compact support (cf. Theorem 4.3.3 
in [ID]), we have 

B\(0) D ch supp(^4/i) = c/tsupp(J-" - (a)) + c/isupp(/i) 

= chsupp(F-\a)) + Bi(0) , 

since Af\ = J- -1 (a) * x /i- Only if < 7 r-1 (a) has singular support {0}, we 
do not obtain a contradiction. If supp(J r_1 (d)) = {0}, then a(z) is a 
polynomial in iz — izi . . . zn (cf. Theorem 2.3.4 in [ID])- 

a) If a(z) =&Gf, then f\{z) = exp{— b} and thus supp(/i) = {0} ^ 

fii(o) (cf. (USD). 

b) If a(z) = ibz with b G M\{0}, then f\{z) = exp{— ibz} and thus 
supp(/i) = {(&,...,&)} 7^ i?i(0). (For 6 = we have case i) a).) 

c) In Theorem|6jin the Appendix, we show that if a(z) is a polynomial 
in iz with exponent larger than 1, then supp(/i) cannot have 
compact support. 

For each of the cases a)-c) we obtain a contradiction which proves that 
the assumption supp(^4/i) C Bi(0) cannot hold for case i). 

ii) That supp(^4/i) C £>i(0) cannot hold if a does not have compact sup- 
port, follows from an estimation of the minimum modulus of a. The 
details are carried out in Theorem |8] in the Appendix. 

In summary, we have shown that supp(^4/i) C B\(0) cannot hold for both 
cases. □ 

Now we are ready to prove the noncausality of G. More precisely 

Theorem 3. Let N G N and a G <S(R ) be a function that is real valued 
and rotational symmetric. Then G defined by nD\) does not satisfy causality 
condition |3|). 

Proof. We perfom a proof by contradiction. Assume that 
(AS) G satisfies condition (J3]) 



8 



and without loss of generality, we can assume that ([Til holds, since J r_1 {a}(x) 
is rotational symmetric. (This can be accomplished by rescaling of space and 
time.) Since G defined by (TfOj) solves evolution equation ( Tl2|) . it follows 

— (•,*) Q chswpp(G(-,t)) C B f (0) 

for each t > 0. In particular, 

(20) cfcsupp(.AG(-, 1)) C Bi(0) . 

But according to Theorem El property (120]) cannot hold. This contradiction 
shows that the assumption (AS) is not true and concludes the proof. □ 



3 A causal model for diffusion 

As shown in Theorem [3l it is not possible that causal diffusion satisfies the 
semigroup property for every time point t > 0. In this section we present a 
model for diffusion that satisfies 

• the semigroup property for a discrete set of time points {r n | n G No}, 

• an inhomogeneous wave equation with a time dependent coefficient 
such that 

• the total mass is conserved and 

• causality holds. 



In the following we use the notion 
Al) c, r G (0, oo) and r m := mr for m G N . 

A2) If t G (r m , T m +i] ( m £ No), then n(t) := m and R(t) :— c(t — T n {t))- 
A3) / : [0, c*] -»■ [0, oo) is such that j[ /(c) dc = 1. 

We emphasize that if t = Tfe, then n(t) = A; — 1 and = cor. Moreover, 
da(x') denotes the Lebesgue surface measure on and |Sr(0)| denotes the 
surface area of the sphere S R (0) C R N . We note that 

J /(x')d ( r(x')= J f(x + Ry)R N - l da(y) and 

( 21 ) 5 fl (x) 5i(0) 

1^(0)1 = 1^(0)1 i^" 1 . 
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3.1 Definition of the model 

First we consider the case where each particular process takes place with the 
same speed c, i.e. c and r are fixed in time. 

Definition 2. Let c, t, r m , n(t) and R(t) be as in Al)-A2). The concentra- 
tion v of a substance diffusing via a constant speed c with initial concentration 
u e L 1 (R N ) is defined by 



(22) 



f v (x' ; r n(t) ) 
v{x,t) = / — — — dcr(x) with v\ t =o = u . 

J \bR(t)\U)\ 



>fl(t)( x ) 



In analogy to wave mechanics, we call A := cr and k := (2tt)/X the wave 
length and the wave number of the diffusion process. 

This definition means that the concentration on the sphere Sfl(t)(x) at 
time t u m determines the concentration in point x at time t. If we consider 
the paths of the particles as lines during the time period (r n (t), T n (t)+i), then 
the particle speed is 

R(t) 



t — T n (t) 



C . 



A more detailed interpretation of this model is given in Subsection 13.31 

Remark 2. For the one dimensional case N = 1, we have Sru)( x ) = {x — 
R(t),x + R(t)}, \S R(t) (x)\ =2 and 

J da(x') = J{5{x'-R(t)) + 5(x' + R(t)))dx' . 

S m (x) H 

Therefore Definition reads as follows 

v ( x , l ) = \ ( v ( x ~ ^(*)' r "(*)) + v ( x + ^(*)' r «(*))) • 

For this special case x i— )■ v(x, t) for each t > is a discrete measure on R. 

That mass is conserved by the diffusion model f)22p follows from the fol- 
lowing lemma. 

Lemma 2. Let u G L 1 (M Ar ) and v be defined as in |Hj]. Then v is well- 
defined and satisfies 

J u(x, t) dx = J w(x) dx for t > . 
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Proof. We perform a proof by induction. 

Let n = 0, i.e. t G (0,r] = (T ,rJ. Then and (J2TJ) imply 

(23) = / ^^ )y) MY) (t6(0,r]) 

Si(0) 

which shows that t>(-,£) is well defined for t G (0, t] and satisfies 
y v(x, t) dx = j J u(x + R(t) y) dx da(y) 

= 1 y w(x) dx . 

Now let nGN, i.e. t G (r n ,r n+ i]. Analogously, (122]) and (1211) imply 

(24) U CM) = y j^^yj d(T (y) (t G (r n ,r n+1 ]) 

Si(0) 

which shows by induction that v(-,t) is well defined for t G (r„,r n+ i]. More- 
over, from the last identity, it follows that 

y u(x, t) dx = y y v(x + y , r n ) dx -i— - d<x(y) 
1 / f(x, r n )dx 



which together with the induction assumption f RN t>(x, t) dx = J RN w(x) dx 
for t G (0, r n ] imply that 

y u(x, t) dx = y u(x, r n ) dx = y m(x) dx. 

R N R N R N 

□ 

Now we define diffusion consisting of various noninterfering processes that 
take place with different speeds c G (0, c*]. 
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Definition 3. Let c*,t* > 0, c£ (0, c*] and r G [r*,oo) 6e snc/i i/iai cr 
c* r* ; / 6e as in A3) and u G i/^R^) fre swc/i i/iai a solution u c oj§ 

c* 

(25) u(x) = J f(c) u c (x) dc (ct = const) 



exists. Moreover, let v c be defined as in t lzfy) with initial data u c . Diffusion 
of an initial concentration u G L 1 (R N ) with speed distribution f is defined by 



(26) u(x, t) — / /(c) v c (jc, t) dc with v\ t=0 = u . 



We call A := c* r* and k := (2tt)/X the wave length and the wave number 
of the diffusion process. 

That model fl26l) conserves mass follows at once from Lemma [2] and prop- 
erty Jo f( c ) dc=1 - 

3.2 The Green function of diffusion and its properties 

Definition 4. Let p(u;x,t) (x G R^, t > 0) denote the concentration of the 
diffusion processes with initial distribution u G L 1 (R 7V ) ; i.e. p(u; x, 0) = u. 

If 

p(u; x, t) — p(5(x); x, i) * x w for each « G L 1 (R ) , 
i/ien we call G(x.,t) := p(5(x); x, i) i/ie Green function of diffusion. 

Below in Proposition [1] and Remark [3] we show that the Green function 
exists if and only if /(c) = 5(c — c ). But before we can do that we need a 
proposition. 

Theorem 4. Let t, c, n(t) and R(t) be as in Al)-A2) and G be defined 
by |Hj] with u = <5(x). Then G satisfied 

(27) G(-,t) = G(;T n{t) )* x G(;t-r n{t) ) for t>0 

and 

G(-,T k+m ) = G(-,T k ) * x G(-,T m ) for k,meN . 



3 This problem is ill-posed and thus cannot be solved without additional information. 
4 This is a modification of the semigroup property of continuous model processes in 
Stochastic Analysis. 
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Proof. By fl2TJ we have for R > 0, x G and # G L 1 (R iV ): 
/ /(x') ^ -Ix-x'Ddx' 

R N 

oo 

= J J f(x + Ry)6(R -\Ry\)P N - 1 

Si(0) 



/(x + i?y) 5(i2o - |fly|) iT" 1 d<r(y) di? 
J f(x+R y)R»- l da(y)= J /(x') da(x') 



51(0) 5 flQ (x) 

From this and (122|) . it follows that 

| MJV G(x', r n ) - |x - x'|) dx' 



\S R (t)(0)\ 



G(x,t) 

In particular, since G(x',0) = 5(x'), 

(28) G(x,t) = g( g^J ) 3 | l) for tG(0,r] 

Combining the last two results yields 



= (G(-,r„) * x G(-,t-r n ))(x) for t > 0, a; G 

Moreover, we get from this for i = r n+ i that 

(29) G(x,r n+1 ) = (G(.,r n )* x G(.,r))(x) 

and consequently by induction G(-,r k+m ) = G(-,r k ) * x G(-,r m ). As was to 
be shown. 

□ 

Proposition 1. Let v be as in Definition^ and G be defined by |HP with 
u(x) = S(x). Then 

u(x,t) = (G(-,t) * x «)(x) t>0, xeR w , 

i.e. G zs i/ie Green function of diffusion with constant speed c. In particular, 
v satisfied 



v(x,t) = J G(x',r n (i))w(x-x / ,t-r„ w )dx / t > 0, x G 



5 This equation is a modification of the Chapman-Kolmogorov equation and since G(-, t) 
is a positive distribution, /it(x) := G(x, t) dx is a measure on M. N . 
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with initial condition t>(x, 0) = u(x). 

Proof. We prove the first statement of the proposition by induction. 
Let n = 0, i.e. t G (0, r] = (r , T\]. According to (1221) and fl2TT) . we have 

(3o> I n^fr dff(y) (ie(M) ' 

Si(0) 

From this and (123 p . we infer for t G (0, t\. 

, . f u(x + R(t)y) , , . 

t,(x ' ()= y 15,(0)1 dg(y > 

%(0) 

/• a(x + fl(t)y-x') ^ , 
-(x) y ^ d.(y)dx 

R N 5i(0) 

= (n * x G(-,t))(x). 
This proves case n = 0. 

Assume that f(x, i) = ((?(•, t) * x tt)(x) holds for t G (r„_ 1 ,r n ] (induction as- 
sumption). Now let t G (r n , r n+1 ]. Analogously as above, it follows from 
that 

v(x + R(t)y,T n (t)) 



.<*,*) = / | S ; (0 <My) 



Si(O) 

u(x',r n (*)) 



gxJ-jgt^y-xQ , 

— isiwi — (y) 



R N Si(0) 

= (v(-,T„( t )) * x G(-,t- r„( t )))(x) 
which together with the induction assumption and (J27j) imply 

t) = [(?(•, r„) * x u] * x i - r n ) = G(-,t) * x u . 

This proves the first claim of the proposition. 

According to (127)) and the first statement of the proposition, we have 

V(; t) = G(-, t)* x U = G(; t - T n ) * x G(; T n ) * x It = G(-, T n ) * x t - T n ) 

which proves the last claim of the proposition. □ 

Remark 3. Let v c be as in Definition [H and v, f be as in Definition [3J 
Moreover, let G c be defined by |HJ|) with u c = 5(x) and 

c* 

G(;t):= [ f(c)G c (;t)dc. 
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From J C /(c) dc = 1 and G c (x, 0) = 5(x), we infer G(x, 0) = 5(x). However, 
in general ( cf. Definition and Definition ^ 

c* 

i;(x,£) = y /(c) (G c (-,0*x w c )(x) dc 



c* 

^ y /(c) (G c (.,£)* x ii)(x) dc=(G(.,0*x«)(x) 

Q 

and therefore we do not call G Green function if /(c) 7^ <5(c — Co). 

The following two propositions are about the smoothness properties of 
the diffusion process defined as in Definition [2] and Definition [31 respectively. 
In the following, we use the subscript c only if necessary. 

Proposition 2. Let v be defined as in Definition^ with u G L l (Jst N ). 

a) Ifu is continuous, then v (x, •) : R — > (0, 00) is continuous for each x G R^ 

and 

u(x,0) := v(x,0+) = w(x) . 

//it is differentiabl^, £/ten§(x, •) : R ->■ (0, 00) is continuous on (r m , r m +i) 
/or eac/i m G N and 



t=T m + 

/ / / 1 

(31) 



dv 



dt 



J t=T„ 



-w*x-^-(-,Tto-) / or mGN . 



In particular, 



dv dv 

-7^{-,T m ) := —(■,T m +) = Q for m £ N . 

Proof, a) Let it be continuous. From ( I23p . it follows that 

|i;(x, £ + s) - u(x,t)| 

< /■ |n(x + + s) y) - i/(x + R(t) y) | 

Si(0) ' * 

and because i?(£) is continuous for £ G (0, r), it follows that t > v(x,t) is 
continuous for t G (0, r). Let £ = r and s G (— r, 0). Then we have (cf. A2)) 

R(t) = c r and R(t + s) = c (r — |s|) 



3 We can also assume Gateaux differentiabilty. 
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and the previous estimation reads as follows 

|v(x, t + s) — v(x,t)\ 

f |u(x + c (r- |s|)y) -n(x + c ry)| 

- J \m\ (y) 

5i(0) 

which shows that 1 1— > v (x, t) is continuous from the left at t — r. 
Let t = and sG (0, r), then ([22]) implies 

u(x,0+) = lim / — dtr(y) 



s^o, s >o y I Si (0)| 

5i(0) 



' /(X) da(y)= V (x,0) 



|5i(0) 



5i(0) 

i.e. t i — y v(x,t) is continuous from the right at t — 0. Moreover, by 

/ i \ i • / w(x + i?(s)y,r) 
u(x,r+)= lim / r^— f| day 

s-+0, s>0 J \bi{0)\ 

dcr(y) = u(x,r) . 



s x (o) 
u(x, r) 



Si(O) 



|5i(0) 



In summary, we have shown that t H- t>(x, t) is continuous for t G (0, r]. 

Assume that t H- t> (x, t) is continuous for £ G (0, r m ]. 

Let t G (r m ,r m+ i). Analogously as above, it follows by induction from 

\v(jc,t+ s) - t>(x, t)\ 

f \v(x + R(t + s)y,T m ) -v(-k + R(t)y,T m )\ 

- J \m\ My) 

Si(0) 

and (HM!) that t i-> f(x, i) is continuous for t G (T m ,r m+1 ]. This proves the 
claim of part a). 

b): Let u be differentiate. From ([Zlj) , it follows 

x /sifoi^^+^Wy' 7 "^*))] -y dcr (y) 

¥ (x ' t} = ^(o)i c ° ' 

i.e. §f(x, ■) is continuous on (r n ( t ), r n (t) + i) and 
(32) ™(. )Tm+ ) = (m G N) , 
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since R(r m +) = lim s ^ 0+ R(r m + s) = and f Sl r \ y dcr(y) = 0. Moreover, 
Proposition [I] implies that 




Figure 1: Left figure: If N — 2 and t G (r, r 2 ), then supp(G(-, t)) is an 
annulus generated by circles of radius Co (t — r) and center at x' G 5 Cor (0). 
Right figure: If N = 2 and t G (r 2 , r 3 ), then supp(G(-, t)) is a disc generated 
by circles of radius c (t — r 2 ) and center at x' G £> cor2 (0). 

Proposition 3. Le£ t> and v c be defined as in Definition^ with u c G L 1 (M. N ) 
for each c G (0, c*]. Moreover, let x G K and £ > 0. 

If (x, c) i — ^ n c (x) continuous on W N x (0,c*) ; t/ien c h-> n c (x, t) zs 
continuous on (0,c*). 

b) If c^r /(c) and c H- u c (x) are continuous on (0,c*) ; iaen t (->■ ||(x, t) zs 
continuous on (r*,oo). 

Proof, a) Let c G (0, c*) and r := A/c. We have to show that 
(34) lim n c _ e (x, t) = lim v c+e ht,t) 

e->0+ e->0+ 
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for the three cases 

t < t , t — T and t > t . 

We prove only the case t = r. The remaining two cases are proven similarly. 
Let c t := c — e and r t := A/c e . From (123|) with t = t < r e and the continuity 
of (x, c) H- ii c (x), it follows that 

hm, Ce (x,i) = hm J My) 

dcr(y) = u c (x,t) . 



Si(0) 
M c (x + ct) 



SiW 



\Si (0)| 



Let c e := c + e and r e := X/c € . From (I23j) and (I24I) with t = r > r e and the 
continuity of (x, c) H- w c (x), it follows that 

hm w Ct x, t )= hm / — — da y 

e^0+ e->0+ _/ 1^1 

Si(0) 

v f f M Ce (x + c £ (t-T £ )yi + c e r e y 2 ) , , , , 




Si(0) 5i(0) 

m c (x + + cry 2 ) 
l^i(0)| 2 

Si(0) 5i(0) 

which proves (|34|) for t = r. 

b) According to Definition [31 we have 



dcr(yi) dcx(y 2 ) = v c (x, r) = v c (x, t) , 



oo 

d(x, £) = y iu(x, £, r) dr with to(x, t,r) = f (^-j u a/t (x, t) 



For t G (t*, oo), let r := r* and r m := mt for m G N. From (1341) . we get 



m-^sezz- = Mt ) Ka(^)];:;:: = q (^n) 



and thus 



^(M) = £ / -^(x, t, r) dr - £ Kx,t, s)] s =Z 

m=0 m=l 

I'm 

oo Tm r +1 dw 

J2 / "^(x,t,r) dr , 



m=0 
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where to(x, t, r) is differentiable for each r G (r m , r m+1 ). This concludes the 
proof. □ 

In the following we focus on the case f(c) = 8(c — cq) and usually write 
G instead of G CQ . We now specify the support of the Green function G and 
show that causality condition fl3]) is satisfied. 

Proposition 4. Letr, c be as in Al)-A2) and /(c) = 5(c— Cq). Moreover, let 
G denote the Green function of diffusion, i.e. G satisfies ( Q?I| ) with u = 5(x). 

a) Let N G N. Then G satisfies causality condition and 

(35) supp(G(;t)) = S Cot (0) for te[0,r\. 

b) IfN> 1, then 

(36) supp{G(;t)) = £ cot (0)V6 Co(2r _ t) (0) for te(r,2r] 
and 

supp{G{-,t))=B cot (0) for t>2r. 

Proof, a) Identity fl35l) follows at once from fT28|) . To prove the rest of the 
claim, we perform a proof by induction. According to f )35|) we have for 
te[0,r]: 

(37) supp(G(.,*)) Cip). 

Assume that (|3T|) holds for t e [0, r m ] (Induction assumption). From (B7|) 
with t G (r m , T m+ i] and fl35l) together with the support theorem of distri- 
butions with compact support (cf. Theorem 4.3.3 in [10J ) , it follows for 

t G (r m , r m+ i] that 

supp(G(-,t)) = c/isupp(G(-,r m )) + ch supp(G(-, t - r m )) 
C B co Tm (0) + B C0 {t _ Tm) (0) = MO) , 

i.e. (137|) holds. By induction, we infer that (|37|) holds for t > 0, i.e. causality 
condition Q is satisfied. 

b) To prove (JM]) we note that ([27]) and (J2SD imply for t G (r, 2 r]: 



G(x,t) = 



/rat ^(c r - 




x'|)5(c (t-r)-\ 


x-x'|)dx' 




^.(0)11^(^(0)1 



Thus, for this case, the support of G(-,t) is the union of the set of spheres 
with radius Co (t — r) and center x' G ^.(O) (cf. Fig. [T] for AT = 2), which 
is the shell £ Cot (0)\£ co{2r _ t) (0). 
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If t G (t 2 ,t 3 ], then according to ( 127)) and f!5oT) supp(G(-, £)) is the union 
of the set of spheres with radius c (t — r 2 ) and center x' G B CoT2 {Q) (cf. 
Fig. [T] for AT = 2), which is the closed ball S Cot (0). By induction, it follows 
from (f271) for £ G (r m , r m+ i] with m > 2 that supp(G(-, t)) is the union of the 
set of spheres with radius cq (t — r m ) and center x' G £> CoT - m (0), which is the 
closed ball £> co t(0). This concludes the proof. □ 

3.3 Interpretation of the diffusion model 

We now give two interpretations of diffusion model (122]) . Firstly, we interpret 
it as transport process and secondly we interpret it as stochastic process. 

Interpretation as transport process 

Let G denote the Green function of the diffusion process f[2"2"j) that takes places 
with constant speed cq. Then the concentration of the diffusing substance is 
given by 

u(x, t) = (G(; t) * x u)(x) u(x, 0) = u(x) 

and thus it is sufficient to explain the process for the case w(x) = 5(x). Let 
n(x) = <5(x). Then the concentration in point x = splits up and starts 
to spread out in each direction at time t — and propagates with speed c 
during the time period (0, r). Because of Lemma [2], mass is conserved during 
this transport process. According to relation (1271) . we have 

G(x,t) = y'G(x,r)G(x-x / ,t-r)dx' for £g(t,t 2 ), 

and hence the concentration on the sphere S CQT (0) = supp (G(x, r)) (cf. ([28]) ) 
splits up and starts to spread out in each direction at time t = t. The con- 
centration propagates with speed c during the time period (r, r 2 ) (cf. Fig. [TJ 
and Proposition H]) such that the total mass is conserved. This process of 
splitting up and spreading out carries on forever if the medium is unbounded. 

We note that this model allows particles to cross (without interfering) and 
consequently not each space-time point can be associated a single velocity 
vector. Indeed, from (J7|) we see that the flux density j(x, r m+1 ) contains the 
velocity vectors 

c = c y for y 6 Si(0) if i)(x + cory,r m )^0. 
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Interpretation as stochastic process 

We can give a probabilistic interpretation of the diffusion model, too. We 
interpret v(-,t) as the probability density of a single particle at time t which 
is part of a diffusion process satisfying (1221) . Mathematically, this means that 

P t {Q) := J u(x',t)dx' 
h 

describes the probability that the particle lies in the region fl at time t. For 
simplicity, we assume N = 3. Let t G (r m , T m +i] and 

1) A t denote the event that this particle arrives (within a fixed cube Q x 
with side length dl) around point x at time t and 

2) B m denote the event that this particle arrives (within a fixed spherical 
shell r Xi?n with thickness di?)0 around surface 5_R( t )(x) at time r TO . 

If event A t occurs, then the probability of the particle to be outside the 
spherical shell r xm at time r m is zero by (|22p . Moreover, if event A t occurs 
and the particle was around point x' G S^(t)(x) at time r m , then the particle 
was around 




for time s G (r m ,r m+ i) . 



However, we do not exactly know where the particle was during this period. 
In Proposition H] we showed that this type of probabilistic process guarantees 
causality and thus in contrast to the noncausal models of stochastic processes, 
the particle lies within a bounded region for each time point. 



3.4 The equation of diffusion 

Although continuity equation (JBJ) with flux density ((7|) can be considered as 
the equation of diffusion, we derive an alternative equation that is a proper 
partial differential equation. We refer only to this equation as the equation 
of diffusion. 

According to Proposition [1] and Theorem HI the concentration of a sub- 
stance diffusing with constant speed, say c , is given by 

v(-,t) = U* X G(-, T n{t) ) * x G(-,t- T n{t) ) t > , 

7 The length dl of the cube determines the thickness dR of the shell by dl 3 = |5di?(x)| di? 
(origin is at x). 
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where u and G denote the initial concentration and the Green function, re- 
spectively. In order to derive the equation of diffusion, we utilize the following 
identities for t G (r m , r m+ i] with m G Nq: 



V 2 <, t) = U * x G(; T m ) * x V 2 G(-, t - T n 



(3f 



(■,t) = U * x G(-,T m ) *, 



d k G 
dt ki 



where k G {1, 2}. We see if (j|rjv x (o iT ) satisfies a partial differential equation 
with time dependend coefficients a(t), b(t),..., then v satisfies the same 
equation with time dependend coefficients a(t — r n ( t )), 6(t — 7Vi(t))> • • • as long 
as t 7^ Tn(t). In addition, we have to pay regard that % has at t = T m a jump 
of size (cf. Proposition |2J) 



(39) 



t=r„ 



<9G 



For the derivation of the equation of diffusion we need two lemmata. 
Lemma 3. Let R(t) be as in A2), then 



5\R(t) 




x|) 




X 





R(t) 



R(ty 



Proof. Without loss of generality we assume that Co = 1, i.e. R(t) = t. Let 
1 ) denote the space of C°°— function with compact support in M JV+1 . 



For i\) G C^{R N+1 ), we have to show that 



J (R(t) - \x\)6'(R(t) - |x|)dxdt 



4(x,t)j^6(R(t)-\x\)dxdt. 



But this is equal to 



^(x,t) (i?(t)-|x|) +V(x,t) 



dx 



fl(t)=|x| 



^(x, |x|) dx 



Riv+i 



which holds for every if; G C£°(R 



□ 
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Lemma 4. Let r, c, r m , R(t) be as in Al)-A2) and let G be defined in §2l 
with u = 5(x). Then g(-,t) := G(-,t — r m ) (m G N ) solves the following 
equation on M N x (r m ,r m+ i) 

1 d 2 g 1 (N-l) dg „ 2 
(40) -A - — -Vo = 

1 J C§ dt 2 + C i?(t) dt 9 



with initial conditions 



(41) g>(x,r m ) = 5(x) and — (x, r m ) = xG 



b) J I^Q) 1UI A C JT! 

From (12"5|) . (j2"Tl) and Lemma [31 we get 



Proof. We note that g(x, t) = G(x, R(t)/co) for x G and t G (r m ,r m+1 ) 



, ( i - |x|) N - 1 5'(i?(t) - |x|) 



V^(x,t) 



\S m (0)\ |x| \S R{t) (0)\ 

5"(R(t) - |x|) 



\S R{t) (0)\ 
Moreover, ([2HD and ([21]) imply 



S'(R(t) - 


-|x|) 


|5jk*)(0) 


|J2(t) 



-(N-l) 



dg J^m-M^ op (N-l) 

and 

&g j 5"(R(t)-\x\) 2<*(N-l)6'(R(t)-\x\) <*(N-1)N 

\S R(t) (0)\ R(t) \S R(t) (0)\ + R(tf 9[ * 

From these identities, we obtain equation f l4U|) for t G (r m , r m+1 ) . 
Since G(x, 0) = 5(x) by (1501). we get 

#(x,r m ) = G(x,0) = 5(x) , 

which proves the first initial conditions in ( 14T|) . The second initial conditions 
in fT4T|) follow from Proposition [51 □ 

To complete the derivation of the diffusion equation we have to incorpo- 
rate the jump condition (13 9 p via a source term (right hand side term). 

Proposition 5. Let cq > be constant, R(t) be as in A2) and v be defined 
as in [22^1 with differentiate u G L 1 (M iV ). Then v satisfies equation^ 



(42) \ 8 2 t v + - . X) d t v - V 2 v = f[u] onR N x(0,oo) 
c cq R(t) 



8 To be precise, we use the generalized time derivative, denoted by dt, since v is not 
differentiable at t — r m . 
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with source term 

1 QQ 

(43) fbA(;t) = ~^2 M *x Yl -^-(^ T -m-)S(t-r n 

and initial conditions 



"° men 



dv 

(44) u(-,0)=ti and — (-,0) = 0. 

L/C 

Proof. All except the source term / follows at once from Lemma |4] with (l38]h 
According to ( 139|) (cf. Proposition E}, fr has a jump at t = r m and thus 



<9t> 



dt 

t=T m + 



dt , 

M *, 



5(t - r m ) 



dt 

5^ ^(-,r m -)5(t-r m ) 



at 2 x ^ at 

where d t and denote the generalized and the pointwise time derivative, 
respectively. From this and (HQ}, we infer equation (142]) with the source 
term (H. □ 



Remark 4. T/ie classical diffusion equation 

— - D V 2 v = with v(x, 0) = u(x) 

can fre obtained via a special limit from model |Hj] /or the time points t = r m 
for m G No with r — > under the side condition 

c 2 t 

D Q := — — = const . 
2N 

Note that D = D(t/2) = \ £ D(s)ds with D(s) defined as in (TJj. For 
these time points, representation [2^1 can be written as follows 

J v(x,r m+1 ) dcr(y) = ^ y [v(x + R{t) y, r m ) + u(x - J?(r) y, r m )] da(y) 

5i(0) 51(0) 

which is equivalent to 

f(x, r m+ i) - u(x, r m ) 

da(y) 



T 

Si(0) 



Am f v(x. + R(T)y,T m ) -2v(x,r m ) + w(x - i?(r) y, r m ) 

iVA) y ^ do(y). 

Si(0) 
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// we perform the limit r — V such that D = const, (and consequently 
c — > oo), then we obtain from 



dv ND f gMx + gy^) * 0Mx,t) 



Si(O) 

Dkj = D 5kj for k, I G {1, 2, . . . , iV}. Apart from the causality prob 

dv 

at 



lem, this limit is not reasonable, since % is not continuous at t = r m (cf. 



Proposition \ 



4 Conclusions 

Since it is not possible that a causal diffusion process satisfies the (strongly 
continuous) semigroup property, we developed a causal model of diffusion 
with constant speed c which satisfies the semigroup property only at a dis- 
crete set M c of time points. In contrast to the classical and noncausal diffu- 
sion model, the causal diffusion model is not differentiable at the discrete set 
M c of time points. This property of non-smoothness forbids the limit pro- 
cess that transforms the discrete and causal diffusion model from Physics and 
Stochastics into the continuous and noncausal standard diffusion model. It 
is surprising that a smoothness assumption, which is too strong, has an effect 
on the causality of a model. Another consequence of causality is that diffu- 
sion (with constant speed) satisfies an inhomogeneous wave equation with a 
time dependent coefficient. We have seen that a diffusion process with vari- 
ing speed c G (0, c*] is a superposition of processes each of which satisfying 
the semigroup property on a discrete set M c such that U c e(0c*] M c = [t*, oo). 
However, the total process does not satisfy the strongly continuous semigroup 
property on [r*, oo) and is not C°° with respect to time. 
Finally, I would like to give two remarks. 

1) It is a valid option to model diffusion- like or stochastic processes by vector- 
valued semigroups. An interesting example from physics may be the Dirac 
equation (cf. [23 EE])- Whether such an approach is more convenient or 
powerful, then modeling with hyperbolic equations or "other equations" is 
not clear. 

2) From the mathematical point of view, it is interesting whether a causal 
variant of the Schrodinger equation is possible such that for example quan- 
tum tunneling obeys causality. It seems natural to the author that the un- 
resolved causality problem of diffusion had an influence on the development 
of quantum mechanics and consequently some of the causality problems in 
this theory may be resolvable. Moreover, it is intriguing that the existence of 
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antiparticles in quantum field theory is based on a "special form" of causality 
(cf. Chapter 2 Section 13 in [21]). 

5 Appendix 

In this appendix we clarify our notion of Fourier transform, convolution, 
formulate the Paley- Wiener Schwartz Theorem which is used for the causality 
analysis, and complete the proof of Theorem [2] from Section [21 

5.1 Definition of the Fourier transform 

The Fourier transform of / £ L 1 {R N ) is defined by 

7-{/}(k) := (27r)-^ 2 J e ik x /(x)dx keR N , 

R N 

where k • x := Ylf=i kj x r I n this notion the convolution theorem reads as 
follows 

nnng) = {2^- Ni2 ns^g} /, g G l\r n ), 

where the space-convolution is defined by 

(/*x0)(x):= f /(x')ff(x-x') dx'. 

5.2 The Paley- Wiener Schwartz Theorem 

We note that the supporting function of a compact set K G is defined by 

N 

(45) iZ*(0 := sup for (Gl". 

Moreover, 

Im(z) := (Im(zi), . . . , Im(zjv)) and \z\ 2 := 2 2 z G . 

i=i 

By -Br(O) we denote the ball of radius i? and center 0. The Paley- Wiener 
Schwartz Theorem (cf. Theorem 7.3.1 in [10]) reads as follows 
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Theorem 5. Let N G N. / : WL N — > R is a distribution with support in the 
compact set K C M. N if and only if 

PI) k G R w i-> J r {/}(k) G R can 6e extended to an entire function denoted 
by zeC N ^ J r {/}(^) G C and 

P^j iaere exzsi M > and C > snc/i t/iat 

TOK^I <C(l + |^|) M e^ (/m(2)) /oroilzeC". 

Tae extension J r {f}(z) is the Fourier- Laplace transform of f, and M can be 
chosen as the order of the distribution f . 

Remark 5. a) If K = [— cq, cq] n , then 



N 

H K (Im{z)) = c \ Im(zj)\ for z G C 



N 



and if K = B Co (0), then 

H K (Im(z)) = c 



N 



S ^Im{zj) 2 for zG 

\ 3=1 



b) LetKC B co (0) and 

(46) w(z 1 ):=(z 1 ,0,...,0)eC n . 
Then from property P2) in Theorem^ it follows that 

logmax| 2l | =r \J r {f}{w{z 1 ))\ 
hm sup ! — ■ < c , 

1 — >oo T 

i.e. F{f} ow is a function of exponential type < Cq (cf. Definition 2.1.3 and 
Theorem 6.8.1 in \T$). 

c) Let /(x) be rotational symmetric in x G R w and J~{f}(z) be entire. Then 
J-{f}{z) is rotational symmetric in z G C^, too, and property P2) of Theo- 
rem^ holds if and only if the following property holds: 

P2') there exist M > and C > such that 

l-HftMzi))! < C(l + \ Zl \) M e HK{Im{zi)) for all z\ G C 

i/ere w is defined as in b). 
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5.3 Final parts of the proof of Theorem [2] 

Let a(z) be entire, G be denned as in (Till and 



h(z) :=G(-,l) = exp{-a(z)}. 



To complete the proof of Theorem [2] in Section [2j we have to prove two 
statements: 

A) If a is a polynomial in iz of degree > 1 (case i) c)), then fi does not 
have compact support, and 

B) if a is entire and J r ~ 1 {a} has not compact support (case ii)), then 
supp(*4/i) C -Bi(O) cannot hold. 

A) a is a polynomial in i z of degree > 1 

In the following a(z) is a polynomial of degree n > 1 (cf. (fl8l) ). Below we 
show that the leading term of the polynomial a(z) is the most important 
part for the problem in question and thus we first investigate the two cases 
of monome with even and odd exponent, respectively. 

Lemma 5. Let N G N, b G R, n G N w and z e C N . If a(z) := bz 2n , then 
fi cannot have compact support. 

Proof. We assume that f\ has support in 



and proof a contradiction. Let Zj = Rj e 1{f: > with Rj > and ipj > for 
1 ^ J < A. For convenience we set 



(47) 



K:= [-L, L] 



N 



for some 



L > 



A' 



(48) 




and 




From a(z) = bz 2n , it follows 



log|/i(*)| 



bRe(z 2n ) 



bRe(R 2n e i2nip ) = -bR 2n cos(2n^). 



We have two cases b > and b < 0. 
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i) Let b > 0. Let (^^) m eN — ((-n > • • • > z n ))meN be the sequence such 
that for j'G{l,..., JV} 



2n (m + tH 



JmeN = (^e 1Vm ) meN with <p 

Then 



7T 

2^ ' 



log|/i(« (m) )| = 4m 2n cos f ^L^j > o 
and (cf. Remark O a) ) 

ifx(Im(z < - m ^)) = L Nm | sin(<^ m )| with lim ip r 

m—^oo 

where sin(y9 m ) ^ 0. Therefore we have for sufficiently large m 
(49) ,„ g | /l(2 «)|« Cl («l^ 2 " 

where Ci, C*2 > are constants. We see that condition P2) in Theo- 
rem is not satisfied for K defined as in fH7|) and consequently, fi does 
not have compact support for b > 0. 

ii) Let b < 0. For the previous sequence (z' m ') m6 N in i) with Lp m := ip : = 
^, we get 

log|/i(>))| = |6|m 2B cos Q) >0 

which leads again to the estimation ( H9|) . Hence condition P2) in Theo- 
rem [5] is not satisfied and consequently fi cannot have compact support. 
This concludes the proof. 

□ 

Lemma 6. Let N G N, b G R, n G and z G C N . If a(z) := ibz 2n+1 , 
then fi cannot have compact support. 

Proof. We assume that f\ has support in K defined as in (J57J) and proof a 
contradiction. Let Zj = Rj e lVj with Rj > and <pj > 0, and R, tp be defined 
as in (HHl) . Similarly as in the proof of the previous lemma, it follows that 

log|/i(*)| = bR 2n+1 sin((2n+ l)tp). 

Again, we have two cases b > and b < 0. 
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i) Let b > 0. Let (z^) me ^ = {{z^, . . . , z ™ )) me rj be a sequence such 
that for j E {1, . . . , N} 



(4 m) ) m eN = (me 1Vm ) m with y? m := yj 



7T 



,r„ W /r« rr„- r- 2 (2 71 + 1) ' 

Then 

log |/i(^ (m) )| = 6m 2n+1 > 
and (cf. Remark a) ) 

H K (\m(z {m) )) = LNm\sm(cp)\ , 

where sin(^) 7^ 0. From this we obtain for sufficiently large m 

'// A (Ii.if:"'")) N 2u + l 



(50) logl/i^l^d 



C 2 



where Ci, C2 > are constants. Hence property P2) in Theorem [5] is 
not satisfied for defined as in (jlTjh 

ii) Let 6 < 0. For the previous sequence (^ m ^) TO eN i n with cp : = 2 (o^i) ' 
we get the estimation (1501) for sufficiently large m, since 

b sin((2n + l)<p) > 0, 

and hence property P2) in Theorem [5] is not satisfied. 

Thus, for both cases, condition P2) in Theorem[5]is not satisfied and therefore 
fi cannot have compact support. As was to be shown. □ 

With the two previous lemmata we now can prove our claim. 

Theorem 6. If P(z) (z 6 <C N ) is a polynomial with degree higher than 1 such 
that f\ = J r_1 {exp(— P)} G S(R N ), then f\ cannot have compact support. 

Proof. Let Q n (z) denote the leading term of — P and g n := exp{Q n }- Then 
the degree n of the leading term is larger than 1 and according to Lemma 
if n is even or Lemma if n is odd, there exists a sequence {z^ m ') me ^ such 
that g n (z^) satisfies the estimation f H9|) or fl50|) with / replaced by g n . This 
together withal 

Re((-P - Q„)(z (m) )) « Re(Q n (z (m) )) for sufficiently large m, 

implies that fi(z^) satisfies the estimation (|49j) or fl50l) . And thus condition 
P2) in Theorem [5] is not satisfied, i.e. f\ cannot have compact support. This 
proves the Theorem. □ □ 



9 This line of argumentation cannot be used for case B) below, where P = a is a power 
series. 
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B) a is entire and T has not compact support 

For the proof we need an estimation of the minimum modulus of entire 
functions. Theorem 3.7.4 in [TJ fulfills our requirements. We state it for 
the case of functions of exponential type r. We recall that F : C — > C 
(N = 1) is of exponential type r if 

log M(r) 
hm sup = r , 

r— too T 

where M(r) := max| z | =r (cf. Definition 2.1.3 in [TJ). 

Theorem 7. If the entire function F : C — >■ C is of exponential type r, then 
for every e, r\ > 0, £ > 1 and sufficiently large R > we have 

log|F(z)| > -A(t + e)£R for \z\ < R 

except in a set of circles the sum of whose radii is at most 2rj^R. Here 
A > depends only on r\ and £. 

Let a(z) be entire and F(z) := exp{— a(z)} be a function of exponential 
type t = 1. Moreover, let w(zi) for z\ G C be defined as in (J15]) . Then 
Theorem [7] implies for the setting e G (0, 1), r\ = 1/8, £ = 2, r = 1 that 

Re(a(w(2i))) < 4Ai? for |zi| = 3R/A > 2t]£R. 

We emphasize that the constant A does not depend on R and therefore R 
can be replaced by any larger number R. For the following proof we use the 
reformulation 

I exp{-a(u;(zi))}| = exp{-Re(a(w(2; 1 )))} 

> exp{— Br} for \zi\ = r 

for sufficiently large r > 0. Here -B := 16 A/3 > and r = 3R/4. 

Theorem 8. Let a and A be defined as in Definitional and G be defined as 
in / TlOj) . If a is entire, J r ~ 1 {a} has not compact support and / TiPj) holds, then 
supp(Afi) C -Bi(O) cannot hold. 

Proof. We perform a proof by contradiction and assume that supp(^4/i) C 
£>i(0) holds. Let w(z 1 ) for z\ G C be defined as in (jlBI) . Then the assumption 
together with Theorem O and Remark O b) imply that 

(52) F{Afi} o w : C -> C, zi i-)> a(tw(zi)) exp{-a(w(^))} 
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is an exponential function of type < 1. We prove the contradiction by showing 
that ^{Afi} o w cannot be of exponential type r 6 (0, oo). Since, by 
assumption, a(z) is entire, rotational symmetric and J r_1 {ci} does not have 
compact support, Remark [5] c) and a) imply that for each C > there exists 
anmeN such that 

max |a(w(zi))| > C max (exp{m |Im(2].)|}) = C exp{m 2 } . 

|zi|=m |zi|=m 

Let M{r) denote the maximum modulus function of J r {Af\}(w(zi)) on the 
circle \z\\ = r. From the previous estimation, (152"]) and fl5"Tl) . we get 

M{m) = max \a(w(zi)) exp{— a(w(zi))}\ 

\zi\=m 

> C exp{m (m — B)} 

for sufficiently large m. Without loss of generality we assume C = 1. Then 
we obtain 

log M(m) 

hm sup > hm sup (m — B ) = oo 

m— >oo fit m—¥oo 

and thus we have proven that F{A f}(w(zi)) cannot be an exponential func- 
tion of type r G (0, oo) which concludes the proof. □ 
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